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Abstract 

Peronia J. Fleming, 1822 is an eupulmonate slug genus with a wide distribution in the Indo-Pacific Ocean. Currently, nine 
species are considered as valid. However, molecular data indicate cryptic speciation and more species involved. Here, we 
present results on a new species found in the Persian Gulf, a subtropical region with harsh conditions such as elevated 
salinity and high temperature compared to the Indian Ocean. Peronia persiae sp. nov. is described based on molecular, 
histological, anatomical, micro-computer tomography and scanning electron microscopy data. ABGD, GMYC and bPTP 
analyses based on 16S rDNA and cytochrome oxidase I (COI) sequences of Peronia confirm the delimitation of the new 
species. Moreover, our 14 specimens were carefully compared with available information of other described Peronia 
species. Peronia persiae sp. nov. is distinct in a combination of characters, including differences in the genital (ampulla, 
prostate, penial hooks, penial needle) and digestive systems (lack of pharyngeal wall teeth, tooth shape in radula, intestine 
of type II). 
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Introduction 

The usually shallow Persian Gulf (average depth 50 m), located in the subtropical northwest of the Indian Ocean, is 
exposed to harsh conditions and limited exchange of water through the Strait of Hormuz. Sea surface temperature 
ranges from 24 to 32°C in the Strait of Hormuz and 16 to 32°C in the extreme northwest, with a short cold spring and 
a long warm summer season (Amini Yekta et al. 2013). Seawater salinity ranges from 37 to 38 ppt in the entrance 
of the Gulf up to 38 to 41 ppt in the extreme northwest, and the range of the tides varies according to the locality 
from 1.4 m up to 3.2 m (Evans 1999). This contrasts the typical water conditions of the Indian Ocean, with sea 
surface temperature ranges of 22.2 to 27.7°C (Sawe 2018), and a lower salinity range of 36.5-37.2 ppt (Pous et al. 
2015). Therefore, the Persian Gulf differs from the Indian Ocean regarding these physiochemical parameters and an 
impact on the diversity and ecology especially of the intertidal fauna is very likely. However, almost no biodiversity 
studies exist from this region. Only a few ones documented the occurrence of marine Heterobranchia along the 
south Iranian shorelines (Amini Yekta et al. 2012; Amini Yekta et al. 2013; Attaran et al. 2015, Rezai et al. 2016), 
with a single species of the eupulmonate Onchidiidae recorded: Peronia peronii (Cuvier, 1804), identified either by 
external morphology (Amini Yekta et al. 2012) or by 18S rDNA sequences in comparison to other heterobranchs 
(Attaran et al. 2015). 

Members of the Onchidiidae, a family within the eupulmonate group Systellommatophora, have a widespread 
distribution and are quite common in the middle and upper intertidal zone, either in rocky, muddy or sandy habitats 
(Dayrat 2009; Zhang et al. 2016). Most species are exclusively found in tropical and subtropical areas of the Indo- 
West Pacific (Wu et al. 2010). Onchidiids have no shell and the notum can be smooth or covered with tubercles 
in different sizes and shapes. Several characters are typical for certain genera, others rather unique. Multiple pho¬ 
toreceptors including dorsal eyes on the notum, for example, are characteristic for the genera Onchidium Buchan- 
nan, 1800, Peronina Plate, 1893, Melayonchis Dayrat & Goulding, 2017, Wallaconchis Goulding & Dayrat, 2018 i 
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Alionchis Goulding & Dayrat, 2018, and Marmaronchis Dayrat & Goulding, 2018, but are absent in Onchidella J.E. 
Gray, 1850, and Onchidina Semper, 1882. Branchial gills on the notum are typical for Peronia. In former times, 
they were also reported for Paraperonia Starobogatov, 1976 (now accepted as synonym of Peronia , Dayrat 2009) 
and Scaphis Starobogatov, 1976. 

Onchidiidae have been recognized since 1815 as a natural taxon and its monophyly was confirmed in recent 
studies (Dayrat 2009; Dayrat et al. 2011). However, the relationships of genera and between species remain poorly 
understood. Furthermore, species identification is hampered by the lack of distinct external features. Although the 
genus Peronia has nine currently accepted species (MolluscaBase, 2018) and is one of the oldest described onchi- 
diid genera, only a few additional studies besides species original descriptions are available. Peronia verruculata 
(Cuvier, 1830) was thoroughly investigated by Awati & Karandikar (1948) using histological methods, but unfor¬ 
tunately, they did not mention the locality of the specimens they investigated. The type locality of P. verruculata 
is mentioned with India, Bandra, and the species seems to have a very wide distribution in the Indo-Pacific Ocean 
(GBIF, 2019). Further information for this species was provided by Plate (1893), Labbe (1934), Britton (1984) and 
Hyman (1999). P. peronii was re-investigated only a few times (Labbe 1934; Plate 1893) and additional information 
on the other seven described species is rarely available. The synonymization of various genera with Peronia has 
rather added to the confusion surrounding this clade, for example Quoyella Starobogatov, 1976 and its monotypic 
species Q. indica Labbe, 1934 (now Peronia indica). In this case, the possession of branchial gills led to the syn¬ 
onymization with Peronia, while other morphological characters i.e. two separate male openings and intestinal type 
five would not support the species in the genus Peronia (Table S1). 

In this study, we address the P. verruculata complex by re-analysing all available literature data and analysing 
all available sequences from NCBI together with new sequences from Iranian specimens. We describe a new species 
from the Persian Gulf based on an integrative morphological and molecular approach and show its distinctiveness in 
morphology from all other valid Peronia species, including the close related species P. verruculata. 


Material and methods 

Sampling was undertaken in the intertidal zone during low tide at Bandar Lengeh (26°33'29"N 54°52'50"E) in 
March 2015, and Lavan Island (26°48'20.99"N 53°16'4.80"E) in February 2016 (Figure 1). In total, 14 specimens 
of Peronia sp. were collected in the intertidal zone from the surface of rocks. Some of them were hidden in rock 
crevices but could be traced by the trails produced during feeding. The animals were photographed alive using a 
digital camera Canon SX160IS and measured for length. For details of locality and preservation method see Table 
1. Temperature, pH value and salinity were measured at Lavan Island in 2016. 



FIGURE 1. Collecting sites of Peronia persiae sp. nov. A. Persian Gulf, general view. B. Rocky shore where the specimens 
were collected during low tides. 
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TABLE 1. Specimens used in this study. Abbreviations: BL=Bandar Lenge; LA=Lavan Island. 


Specimen 

Preservation 

Purpose 

Length of preserved animal 
(mm) 

BL 1 

ETOH 96% 

Paratype; DNA barcoding, SEM, dissection 

37 

BL 2 

Formalin 

Histology 

32 

LA 1 

ETOH 96% 

Paratype; DNA barcoding, SEM, dissection 

32 

LA 2 

ETOH 96% 

Paratype; DNA barcoding, SEM, dissection 

22 

LA 3 

ETOH 96% 

Voucher; DNA barcoding, SEM, dissection 

26 

LA 4 

ETOH 96% 

Voucher; DNA barcoding, SEM, dissection 

28 

LA 5 

ETOH 96% 

Voucher; DNA barcoding, SEM, dissection 

34 

LA 6 

ETOH 96% 

Voucher; DNA barcoding, SEM, dissection 

32 

LA 7 

ETOH 96% 

Holotype; DNA barcoding, SEM, dissection of 

head area 

35 

LA 8 

ETOH 96% 

SEM, DNA barcoding, dissection 

25 

LA 9 

ETOH 96% 

SEM, DNA barcoding, dissection 

22 

LA 10 

ETOH 96% 

SEM, DNA barcoding, dissection, Histology of 
genital system 

13 

LA 11 

Formalin 

Histology 

29 

LA 12 

ETOH 96% 

Micro CT, dissection 

31 


Anatomical observation. Dissection was performed under an Olympus SZX12 stereo microscope. Radula, 
penis, and needle of penial accessory gland were extracted by using a 5% KOH solution, dried and analysed using 
a Scanning Electron Microscope (ZEISS Sigma VP 300). For Micro-CT analyses one specimen was stained for two 
days in 1% iodine dissolved in 100% Ethanol (I2E). Reconstructed images were analysed in CTVox. For histologi¬ 
cal analyses, two specimens and the genital system of an additional specimen were dehydrated in EtOH and embed¬ 
ded in hydroxyethyl methacrylate (Heraeus Kulzer GmbH) for serial sectioning. Sections (2.5 pm) were stained 
with toluidine blue, photographed subsequently under a ZEISS Microscope (Imager.Z2m) and analysed with ZEN 
software (ZEISS) at Alexander Koenig Research Museum in Bonn, Germany. 

DNA extraction, PCR, and DNA sequencing. DNA isolation was carried out using the Qiagen DNeasy Blood 
and Tissue kit, following manufacturer’s instructions. Partial sequences of mitochondrial COI (ca. 680 bp) and 
ribosomal 16S (ca. 550 bp) were amplified by polymerase chain reaction (PCR) using the primers LCOI490-JJ 
(5’-CHACWAAYCATAAAGATATYGG-3’) and HC02198-JJ (5’-AWACTTCVGGRTGVCCAAARAATCA-3’) 
(Astrin & Stiiben 2008) for COI; 16Sar-L (5’-CGCCTGTTTATCAAAAACAT-3’) and 16Sbr-H (5’- CCGGTCT- 
GAACTCAGATCACGT-3’) (Palumbi et al. 1991). The thermoprofile during the PCR used for COI and 16S was: 
15 min at 95°C; 40 cycles of reaction conditions (48 cycles for 16S) involved an initial denaturation of 94°C for 
35sec (45 sec for 16S), subsequent annealing at 55°C for 90 sec (56°C for 45 sec for 16S), elongation at 72°C for 
90 sec and final elongation step of 72°C for 10 min. Sequencing was performed by Macrogen Europe (Amsterdam, 
Netherlands). Sequences are deposited in GenBank with the accession numbers listed in Table 2. 


TABLE 2. Name of species, locality information, COI and 16S GenBank accession numbers for genus Peronia. # Regis¬ 
tered but unpublished data of specimens taken from NCB1. 


Preliminary name taken from 
literature (NCBI) 

New assigned / 

Confirmed name 

Locality 

GenBank COI 

GenBank16S 

Scaphis sp. 


Philippine Islands 

HQ660050 

HQ659918 

Onchidium verrucosum 

Peronia verruculata 
(misspelling by Gobbe- 
ler and Klussmann-Kolb 

2011) 

Australia 

(Queensland) 

EF489391 

EF489316 

Peronia sp. 1 


Hawaii 

HQ660038 

HQ659906 
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TABLE 2. (Continued) 


Preliminary name taken from 

New assigned / 

Locality 

GenBank COI 

GenBank16S 

literature (NCBI) 

Confirmed name 




Peronia sp. 2 


Oman 

HQ660044 

HQ659912 

Peivnia sp. 3 


Queensland (Aus¬ 
tralia) 

HQ660048 

HQ659916 

Peronia sp. 4 


Mozambique 

HQ660045 

HQ659913 

Peronia sp. 5 


Mozambique 

HQ660047 

HQ659915 

Peronia sp. 6 


Indonesia (Su¬ 
lawesi) 

HQ660046 

HQ659914 

Peronia sp. 7 


Hainan (China) 

HQ285979 

HQ285967 

Peronia sp. 8 


Hainan (China) 

HQ285980 

HQ285968 

Peronia sp. 9 


Hainan (China) 

HQ285981 

HQ285966 

Peronia sp. (unpublished 

Peronia verruculata 

Hainan (China) # 

JN543165 

JN543101 

paper by Chen et al. 2011) 

(reassigned by Sun et al. 
2014) 




Peronia cf. verruculata 

Peronia sp. (reassigned by 

Okinawa (Japan) 

HQ660043 

HQ659911 

(Dayrat et al. 2011) 

Peronia peronii 

Dayrat et al. 2016) 

Guam 

HQ660041 

HQ659909 

Peronia cf. peronii 


Mozambique 

HQ660042 

HQ659910 

Peronia verruculata 


Fujian (China) # 

GUI 66568 

- 

(unpublished paper by Wang 


Fujian (China) # 

GU166567 

- 

et al. 2009) 


Fujian (China) # 

GUI 66566 

GQ985284 



Fujian(China) # 

GUI 66564 

GQ985282 



Fujian(China) # 

GU166563 

GQ985281 



Hainan(China)# 

GUI 66561 

GQ985277 



Hainan(China) # 

GU166559 

GQ985275 



Hainan(China) # 

GU166558 

GQ985274 



Hainan(China) # 

GUI 66562 

GQ985278 



Hainan(China) # 

GUI 66560 

GQ985276 



Fujian(China) # 

GUI 66565 

GQ985283 



Hainan(China) # 

GU166557 

GQ985273 

Peronia verruculata (unpub- 


Fujian (China) # 

JN543154 

JN543090 

lished paper by Chen et al. 


Hainan (China) # 

JN543153 

JN543089 

2011) 


Zhanjiang (China) # 

JN543152 

JN543088 

Peronia sp. 


Singapore 

MH002607 

- 




MH002590 

- 




MH002605 

- 




MH002603 

- 




MH002600 

- 




MH002599 

- 




MH002594 

- 




MH002591 

- 




MH002589 

- 




MH002586 

- 




MH002585 

- 




MH002580- 

- 




MH002582 

MH002575- 

MH002578 

- 

. continued on the next page 
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TABLE 2. (Continued) 


Preliminary name taken from 
literature (NCBI) 

New assigned / Locality GenBank COI GenBank 16S 

Confirmed name 

Peronia sp. 2 

MH002592 

MH002596 

MK142731 

MK142730 

MK142725- 

MK142728 

MK142720- 

MK142723 

MK142717 

MK142715 

MK142714 

MK142712 

MK142704 

MK142706- 

MK142709 

MH002597 

MH002604 

MH002579 

MH002598 

MH002602 

MH002583 

MK142710 

MH002593 

MK142729 

MK142713 

MH002595 

MK142732 

MK142724 

MK142716 

MH002588 

MK142705 

MK142719 

MK142711 

MH002584 

MH002606 

MH002601 

MK142718 

Singapore MK142680- 

MK142703 

MH002574 

MH002573 

MH002569 

MH002567 

MH002566 

MH002559- 

MH002563 
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TABLE 2. (Continued) 


Preliminary name taken from New assigned / Locality GenBank COI GenBank 16S 

literature (NCBI) Confirmed name 




MH002548- 

- 



MH002556 

MH002546 




MH002545 

- 



MH002571 

- 



MH002572 

- 



MH002547 

- 



MK142694 

- 



MK142689 

- 



MH002564 

- 



MH002550 

- 



MH002557- 

- 



MH002558 

MH002570 


Peronia persiae sp. nov. 

Bandar Lenge 

MK312166 

MK312167 


Lavan Island 

MK993386- 

MK993398- 



MK993395 

MK993407 

Peronia sp. 7 

Indonesia, Bangka 
Island 

MK993396 

MK993397 


Phylogenetic reconstruction. Sequences were edited using BioEdit (ver.7.2.6.1) (Hall 1999) and aligned using 
MAFFT (Katoh et al. 2002) in Geneious v7.1.9 (Kearse et al. 2012). Sequences of other species from all available 
genera of the family Onchidiidae (446 sequences) and outgroup (three sequences) were downloaded from Gen¬ 
Bank, from Wang et al. (2009), Chen et al. (2011), Dayrat et al. (2011, 2014, 2016, 2017, 2018, 2019), Dayrat & 
Goulding (2017); Cumming et al. (2014), Goulding etal. (2018a, 2018b, 2018c) and Layton et al. (2014) (Table 2). 
All sequences were included in the phylogenetic analysis to evaluate the affiliation and monophyly of Peronia for 
subsequent analyses. Two species of Siphonaria G. B. Sowerby I, 1823 (.S', zelandica Quoy & Gaimard, 1833 and 
S. sirius Pilsbry, 1894), considered as the closest relatives to the Onchidiidae (Wagele et al. 2014) were included 
as outgroup for rooting the tree. After trimming, the alignments comprised 454 bp for 16S, 635 bp for COI and 
1075 bp for the concatenated dataset. The maximum likelihood (ML) analysis was run in IQ-TREE (Nguyen et al. 
2014; Trifmopoulos et al. 2016) based on the concatenated data set using the online version 1.6.3 on a Webserver 
(http://iqtree.cibiv.univie.ac.at/), with the GTR model for each gene. Support values were calculated based on 1000 
ultrafast bootstrap replicates and the approximate likelihood ratio test (SH-aLRT) (2000 replicates). Subsequent 
analyses were performed with a reduced and concatenated alignment, only including Peronia sequences (Table 2), 
with identical sequences removed. For this analysis, Wallaconchis graniferus (Semper, 1880) was chosen as out¬ 
group, based on the results of the general Onchidiidae analysis. Species delimitation tests were performed with the 
same taxa, but with 16S and COI genes datasets separated, and the results presented on the COI (Figure 11) and the 
concatenated trees (Figure 12). Dendroscope (version 3.5.8) (Huson & Scornavacca 2017) and Inkscape (version 
0.92) (https://inkscape.org/en/) were used to edit the phylograms. 

Species delimitation. The Automatic Barcode Gap Discovery (ABGD) (Puillandre etal. 2011), General Mixed 
Yule Coalescent model (GMYC) (Monaghan et al. 2009; Pons et al. 2006) and Poisson Tree Processes (PTP) 
(Zhang et al. 2013) were applied for delimiting the species. ABGD is independent of predefined species groups 
(Puillandre et al. 2011; Padula et al. 2014) and was applied to COI and 16S datasets separately and as concatenated 
datasets, using default values under the Kimura, K80model. GMYC is a likelihood method for delimiting species 
by fitting within- and between-species branching models on phylogenetic trees (Pons et al. 2006; Monaghan et al. 
2009). bPTP is an updated version of the original maximum likelihood PTP program (maximum likelihood PTP 
search result is part of the bPTP results) and adds Bayesian support (BS) values to delimited species on the input 
tree (Figures S1-S2) (Zhang et al. 2013). Ultrametric starting trees for GMYC were generated using penalized 
likelihood (Sanderson 2002) using the chronos function in the APE package (Paradis et al. 2004) in R. Four clock 
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models: strict, discrete rate variation with ten rate categories, correlated and uncorrelated-relaxed were fitted on the 
ML trees for COI and 16S using the chronos function in the APE package (Paradis et al. 2004). The best model was 
selected using the phi information criterion, which takes the penalized term into account (Paradis 2013). All models 
were fitted for trees based on both COI and 16S with the smoothing parameter, lambda, set to 0.1, 0.5 and 1.0, and 
in each case, the strict clock was found to be the best model. Ultrametric trees for both genes based on a strict clock 
with lambda=l were used in subsequent analyses. Single threshold GMYC analyses were performed in R using the 
SPLITS package (Ezard et al. 2009). The PTP model assumes that the number of substitutions is significantly higher 
between species than within species, which is reflected in the branch lengths (Zhang et al. 2013). bPTP analyses 
were run on COI and 16S trees resulting from the ML analyses and uploaded, individually, as Nexus files to the 
PTP web server (http://species.h-its.org/) (Zhang et al. 2013). Trees were rooted and W. graniferus was included 
as outgroup. bPTP graphic results for each gene were presented as PhyloMaps (Zhang et al. 2011). COI minimum 
and maximum pairwise uncorrected p-distances, between and within the main clades, were calculated using Species 
Identifier (Meier et al. 2006). 


Results 

Systematics 

Family Onchidiidae Rafinesque, 1815 
Genus Peronia J. Fleming, 1822 

Type species Peronia peronii (Cuvier, 1804) 

Diagnosis. Presence of contractible gills distributed irregularly at least in the posterior part of the dorsal notum. 
Notum always narrower than the pedal sole; warty with dorsal eyes single or in groups. Some species character¬ 
ized by spicules in the integument. Intestine usually of type I or type II. Rectal gland absent. Penis short with many 
spines. Penial gland very long with a large muscular sac. Needle apparatus in penial gland sac present (Labbe 1934; 
Britton 1984; Dayrat2009). 


Peroniapersiae sp. nov. Maniei, Espeland, Movahedi & Wagele 

(Figs. 2-10) 

ZooBankLSID: um:lsid:zoobank.org:pub:2F2B0734-03E2-4D94-A72D-9E43A132DlDE 

Type material. Holotype: Zoologische Staatssammlung Miinchen (ZSM Mol 20180017), Lavan Island 
(26°48'20.99"N 53°16'4.80"E), collected in February 2016, 35 mm in length. Paratypes: Zoologische Staatssam¬ 
mlung Miinchen (ZSM Mol 20180018), Bandar Lengeh (26°33'29"N 54°52'50"E), collected in March 2015 (1 
specimen), 32 mm in length; Lavan Island (26°48'20.99"N 53°16'4.80"E), collected in February 2016 (2 speci¬ 
mens), 22 and 37 mm in length. 

Etymology. Peronia persiae sp. nov. is named after the home country of the first author, Iran (Persia). 
External morphology (Figure 2). Living specimens with 20-65 mm in length, and 13-37 mm when preserved 
in formalin. Elongate while moving, but more oval to round in outline when resting. Notum of living animals muddy 
green to grey in colour and covered with tubercles (papilla) (Figure 2A-B). Foot elongate, smooth and light green; 
covered by hyponotum (Figure 2A). Eight to 16 tubercles bearing dorsal eyes in groups of two to four (Figure 2C). 
In posterior region of the mantle, six to 16 irregularly branched branchial gills (so called gill-trees) (Figure 2D). 
Gills expanding only when the animal is submerged. Two retractable tentacles with an apically lying eye between 
mantle and foot in the head region. Ventrally lying mouth surrounded by oral lips. Male genital pore between right 
tentacle and labial palp (Figure 2A). Anus located ventrally at the posterior end of the body in the middle between 
mantle and foot; female opening on the right side of the anus; pneumostome posterior to the anus and opening only 
when the slug is out of the water. 
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FIGURE 2. Peronia persiae sp. nov. External characters. A-B. Ventral and dorsal view of Peroniapersiae sp. nov. C. Tuber¬ 
cles with dorsal eyes. D. Branchial gill at the rear of the back. Abbreviations: A=anus; Bg=branchial gill; De=dorsal eye; F= 
foot; Fa=female genital aperture; Fln=hyponotum; Lp=labial palp; N=notum; Pn=pulmonary aperture; Se=stalk eye; T=tentacle 
bearing eye; Tb=dorsal tubercle. 

Integument. Epidermis composed of small cells; completely covered by a thin cuticular-like layer. Many bluish 
stained glandular cells present, reaching deep into the notum tissue; glandular cells staining violet (mucopolysac¬ 
charides) interspersed. Some cells without distinct contents (probably empty gland cells), also reaching into notum 
tissue (Figure 3A). 

Vision system. Oval shaped stalk eye on the tip of each tentacle; consisting of cornea, lens and pigmented 
retina layer (Figure 3B). The latter comprising a villous, pigmented layer, a somatic layer and neural layer with 
nerve fibres. The dorsal eyes on the notum of inverse nature: the position of the retinal layer inverted with the main 
nerve penetrating the pigment layer surrounding the retina. Lens composed of two parts: Principal lens consisting 
of one cell with a large nucleus and thick microvilli layer facing towards epidermis. Accessory lens lying beneath 
principle lens; formed by at least one cell with a large nucleus. Epidermis covering dorsal eyes similar to surround¬ 
ing epidermis, but with smaller cells above the eyes (Figure 3C). Photoreceptors as extraocular structures scattered 
beneath the epidermis, partly single or forming clusters of up to eight receptors and sometimes in close vicinity to 
dorsal eyes. Each receptor composed of a cell with a thick layer of microvilli and large nucleus. No retina or cornea 
observed. Photoreceptors separated by connective tissues and probably muscle cells (Figure 3D). 
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FIGURE 3. Histological sections of epidermis and sensory organs of Peronia persiae sp. nov. A. Cross section of the mantle. B. 
Cross section of stalk eye. C. Longitudinal section of dorsal eye. D. Cross section of four dermal photoreceptors. Abbreviations: 
Al=accessory lens; E=empty unicellular gland; EP=epidermis; G=glandular vesicle consist of mucopolysaccharides (stained 
violet); U=unicellular gland; Mv=microvillus; Nu=nuclei; Op=optic nerve; Pg=pigmented layer; Pl=principal lens; Ps=pseudo 
cornea; Rl=retina layer. 

Digestive system. Radula broad, with teeth lying on a thick cuticular membrane (Figure 4A). Interior of odon- 
tophore filled with red staining substance reminding of hyaline material or connective tissue (“elastic non-cellular 
substance”, Awati & Karandikar 1948, p. 14). Odontophore flanked by the lateral cartilage-like radula cushions as 
well as pharyngeal muscles; radula cushions connected by a rigid area probably consisting of red stained connec¬ 
tive tissue (Figure 4A). Radular formulae counted for 11 specimens; formulae ranging from 49 x 47.1.47 (specimen 
length alive 22mm) up to 71 x 87.1.87 (specimen length alive 65 mm). Rachidian teeth tricuspid with a main middle 
cusp and one distinctly narrow and long lateral cusp on both sides of main cusp (Figure 5A). Height of rachidian 
teeth about 50 pm. First inner lateral teeth smaller than other laterals (Figure 5B); height of unicuspid lateral teeth 
gradually increasing from about 60 pm to 75 pm and decreasing again towards the lateral rim. Lateral teeth with a 
conical shape seen from lateral view and inner side curved in a concave way (Figure 5C). Tips (seen from above) 
usually broadened and with a blunt, spatulate apex (Figure 5A-B). Salivary glands on both sides of the oesophagus 
with small ducts leading into pharynx close to the transition into the oesophagus; composed of many finger-like 
tubes combining into clusters and fonning a grape-like structure. Glandular cells with reddish to violet staining 
granular contents (Figure 4B). Oesophagus with highly folded epithelium composed of ciliated columnar cells, 
covered in many areas with a thin homogenously staining (probably cuticular) layer—sometimes more greenish, 
sometimes bluish (Figure 4C). Oesophageal folds underlain with red stained connective tissue and surrounded first 
by a longitudinal and then by a circular muscular layer. Oesophagus entering the first part of the stomach. Stomach 
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consisting of four parts. First part characterized by a thin epithelium; receiving the ducts of the dorsal and left lateral 
lobes of the digestive gland. Second chamber strongly muscular (Figure 4D), swollen and stratified; receiving the 
duct of the posterior lobe of the digestive gland. Third chamber funnel-shaped; pigmented on the outer side; its epi¬ 
thelium forming a highly folded structure with dendritic branches filling nearly completely the internal lumen (Fig¬ 
ure 4E). The last chamber representing a small widened section at the beginning of intestine with only thin ridges 
internally and ciliated cells in the epithelium. Intestine forming two distinct long loops lying close together (Figure 
6A—C), thus fitting best type II according to the definition of Labbe (1934). Epithelium with light blue stained cells 
and with violet stained goblet cells (Figure 4F). Digestive gland composed of many lobes; epithelium exhibiting 
cells which excrete substances stained in various shades of violet to blue (Figure 4G). 



FIGURE 4. Histological sections of the digestive tract of Peronia persiae sp. nov. A. Cross section through the posterior region 
of the pharyngeal cavity and radula. B. Salivary gland cells. C. Cross section of oesophagus. D. Second chamber of stomach. 
E. Third chamber of stomach with highly folded (branched) epithelium. F. Cross section of intestine. G. Cross section of diges¬ 
tive gland. Abbreviations: Ce=columnar epithelium; Cl=cuticular lining; Ct=connective tissue ; Ep=epithelium; Et=non-cel- 
lular elastic tissue; Gc=goblet cells; Pc=pharyngeal cavity; R=radula bearing teeth; Rs=radula support; Tm=transverse muscle; 
Lm=longitudinal muscle. 
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FIGURE 5. Scanning electron micrographs of the radula of Peronia persiae sp. nov.: A. Rhachidian teeth. B. Lateral teeth. C. 
Lateral tooth. 




FIGURE 6. Digestive system of Peronia persiae sp. nov. A. Dorsal view. B. Dorsal view with digestive gland removed. C. 
Ventral view. Abbreviations: ddg=dorsal lobe of digestive gland; i=intestine; oddg=opening of dorsal lobe of digestive gland; 
oldg=opening of lateral lobe of digestive gland; opdg=opening of posterior lobe of digestive gland; pdg=posterior lobe of diges¬ 
tive gland; r=rectum; stl=stomach chamber 1; st2=stomach chamber 2; st3=stomach chamber 3; st4=stomach chamber 4. 
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FIGURE 7. Histological sections of the genital system of Peronia persiae sp. nov. A. Ampulla. B. Cross section of penis. C. 
Penial accessory gland duct. D. Female gland mass (mucus gland and albumen gland). E. Receptaculum seminis. F. Spermath- 
eca. Symbols: arrow=penial hook; asterisk=eggs; dot=sperms; head arrow=muscular layer. 


Reproductive system. Hermaphrodite gland (gonad) compact, whitish with bright dots representing the oogo- 
nia; located slightly on the left side of the visceral cavity, next to the diaphragm. Hermaphroditic duct originating 
from gonad and continuing to a widened area, possibly the ampulla, containing sperm, as well as a few oocytes; duct 
with small epithelial cells (Figure 7A, Figure 9A). Vas deferens fonning a long tube starting from hermaphroditic 
duct in the posterior third of body, running to the right side of the head, and ending in the penis; surrounded with 
a thick muscle layer and lined internally by ciliated cells. Prostate gland not observed. The anterior organs of male 
system consisting of the penial structure (penis with vas deferens and attached retractor muscle ending in penial sac) 
and penial accessory gland (glandular duct, hollow needle and end sac of penial gland). Both complexes opening 
into a common vestibule and sharing same male opening (Figure 9B). Penis usually inverted into the vas deferens 
(Figure 7B) ending in the penial sac. The outer epidermis of penis covered with spines. During evagination, the 
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usually outer layer turns to the outer side (Figure 8A). The hooks at the apical and distal areas are short and coni¬ 
cal while the rest of spines in the middle part are long, sharp and curved downward at the tips. Fork-shaped spines 
existing in different areas of the penis (Figure 8B). Retractor muscle attached to the base of the penis in front part, 
running far into the posterior part of body. Penial accessory gland duct long, heavily coiled and swollen in posterior 
part, close to the opening; epithelium composed of columnar cells with apocrine secretion (Figure 1C). Light yel¬ 
low needle apparatus straight, narrow and hollow; in one case tip slightly curved; length of needles around 1.3 mm 
(Figure 8C-D). The opening of needle usually blunt, only in one individual of sickle shape (Figure 8E). Sac-like 
structure without papillae. Proximal spermoviduct folded and embedded within the female gland mass. The latter 
composed of capsule gland (albumen gland), followed by membrane gland, leading into mucus glands (Figure 7D, 
Figure 9A). Membrane gland located ventrally underneath capsule gland, composed by columnar cells filled with 
reddish granules, mucus gland composed of columnar cells containing basal nucleus and many small, ovoid, violet 
staining granules. Presence of spiral glands could not be verified. Whitish receptaculum seminis composed of small 
epithelial cells, with attached sperm heads (allosperms) (Figure 7E). Brownish round to elongate spermatheca (or 
bursa copulatrix in the sense of Schmekel 1971) connecting to short distal oviduct; greenish staining contents pres¬ 
ent inside spermatheca; epithelium composed of apocrine secreting cells (Figure 7F). Female opening near the anus 
slightly to the right side. 



FIGURE 8. Cuticularized structures in the genital system of Peronia persiae sp. nov. A-C. Scanning electron microscopy. A. 
Partly everted penis. B. Forked penial hook. C. Typical needle in penial structure. D. Light microscopy. Curved needle, observed 
only once. E. Shape of needle aperture retrieved from literature and current study: (a) P. gondwanae (Labbe 1934); (b) P. indica 
(Labbe 1934); (c) P. verruculata (Awati & Karandikar 1948); (d) P. verruculata (Plate 1893); (e) P. branchifera (Plate 1893); P. 
persiae sp. nov. (drawing and SEM picture). 
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FIGURE 9. Reproductive and excretory system of Peronia persiae sp. nov. A-B. Female part of the hermaphroditic genital sys¬ 
tem. B. Male copulatory parts. C. Posterior part of the body opened and displayed from dorsal side. Abbreviations: A=ampulla; 
Au=auricle; Di=diaphragm; Fgm=female gland mass; FIg=hermaphrodite gland; Ne=nephridium; Od=oviduct; P=penis; 
Pb=pulmonary cavity branches; Pgd=penial gland duct; Pgs=penial gland sac containing needle; Ps=penial sac; Psh=penial 
sheath; Rm=retractor muscle; Rs=receptaculum seminis; Sod=spermoviduct; Sp=spennatheca; V=vagina; Vd=vas deferens (to¬ 
wards head); Ve=vestibule; Ven=ventricle. Connection from spermoviduct to female gland mass not clearly seen and indicated 
here with dotted line. 

Circulatory and excretory system. Heart on the right side, divided into an anterior ventricle and posterior au¬ 
ricle within pericardial cavity. Nephridium posterior to the diaphragm. Extension of nephridium larger on the right 
side, than on the left side (Figure 11C). 

Phylogeny. The phylogenetic analyses of the concatenated dataset (Figure 10) covering 11 of the 16 acknowl¬ 
edged genera and including 457 sequences, provide evidence of the monophyly of Peronia with high bootstrap 
support in the ML analyses (99) and close relationship with a clade formed by the monophyletic Wallaconchis 
(96), Alionchis (100) and Paromoionchis Dayrat & Goulding, 2019 (100). The close relationship of the four genera 
is supported by a bootstrap value of 97. One specimen retrieved from NCBI as Scaphis sp. groups within the ge¬ 
nus Peronia. Further monophyletic genera are Peronina (100), Melayonchis (76), Marmaronchis (100), Onchidina 
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(100), and Onchidella (99). Monophyly of Onchidium, as shown in Dayrat et al. (2016) with a reduced dataset of 
Onchidiidae, cannot be seen in our tree, however misidentification of the paraphyletic genus Platevindex H. B. 
Baker, 1938 with members grouping with Onchidium cannot be excluded. In this overall analysis, P. persiae sp. 
nov. is clearly separate from all other Peronia species (Figure 10). These results are confirmed by analysing the 
reduced datasets including only Peronia sequences and running the analyses with COI only (Figure 11, Table 3), or 
the concatenated alignment with COI and 16S, and analyses of distance values (Figure 12). In these two analyses, 
the 11 specimens of the new species always group as a separate clade, usually as sister of Peronia sp. 3, and close 
to Peronia sp. from Singapore (COI, Figure 11) or close to Peronia sp. 4 and sp. 5 from Mozambique (concatenated 
data set, Figure 12). 
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FIGURE 10. Phylogenetic reconstruction of Onchidiidae resulted from Maximum likelihood (ML) analysis of concatenated 
data set (COI and 16S) with Siphonaria as outgroup. Peronia persiae sp. nov. highlighted in bold. Numbers before and after 
slash indicate approximate likelihood ratio test (SH-aLRT) and ultrafast bootstrap values respectively. 
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FIGURE 11. Maximum likelihood tree of the genus Peronia based on COI data set, Wallaconchis graniferus used as outgroup 
(not shown). Groups resulting from ABGD, GMYC, and bPTP tests indicated on the right side. Numbers before and after slash 
indicate approximate likelihood ratio test (SH-aLRT) and ultrafast bootstrap values respectively. 


Species delimitation. All three species delimitation methods applied to the COI (Figure 11) and 16S datasets 
indicate nine well supported groups in the genus Peronia (Figures 11-12, S1-S4). The Iranian specimens always 
form a separate taxon from P. verruculata specimens or any other Peronia species. The result of the ABGD test 
based on COI data shows a single and obvious barcode gap: intraspecific variability is less than 2% and interspe¬ 
cific variability more than 4%. The distance between the new Iranian and other close related species, Peronia sp. 
4 and Peronia sp. 5, ranges between 7.0-8.0% (Table 3). All three tests applied on the separate two genes (Figure 
12) show the same results: one Singapore group, assigned to Peronia sp., is indicated as a separate species. Another 
group, assigned to Peronia sp. 2 from Singapore and Oman, together with sequences assigned to P. verruculata and 
Onchidium verrucosum (misspelling of P. verruculata) from China and Australia, and a few sequences from differ- 
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ent localities (Peronia sp. from China, Peronia sp. 6 and Peronia sp. 7 from Indonesia, Scaphis sp. from Philippine 
Islands) are considered to be a single species. Results of GMYC and bPTP analysis, using the 16S data set, shows 
Peronia sp. 7 from Indonesia and Scaphis sp. as a separate group. Peronia sp. 4 and sp. 5 from Mozambique are 
indicated also as a single and distinct species; however, this last result was not retrieved in the GMYC and bPTP 
analysis using the 16S data set (Figure 12). Peronia sp. 3 from Australia and Peronia sp. 1 from Hawaii also are dis¬ 
tinct species in all analyses. Peronia cf. peronii from Mozambique is considered as a distinct species, when applying 
only the COl dataset. This result is not confirmed in the GMYC and bPTP analysis using only the 16S data set. P. 
peronii from Guam (indicated as a separate species) clusters with a sequence of Peronia cf. verruculata from Japan 
(Peronia sp. according to Dayrat et al. 2016), which also forms a separate species (confirmed in all analyses). 
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FIGURE 12. Maximum likelihood tree of the genus Peronia based on the concatenated alignment (COI and 16S). Groups 
resulting from ABGD, GMYC, and bPTP tests indicated on the right side. Numbers before and after slash indicate approximate 
likelihood ratio test (SH-aLRT) and ultrafast bootstrap values respectively. 
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TABLE 3. Intra- and interspecific pairwise genetic distances. Ranges from minimum to maximum distances are indicated 
(as percentages). * indicating a group of sequences comprising Peronia verruculata (China), Onchidium verrucosum 
(Australia), Peronia sp.2 (Oman), P. sp.6 (Indonesia), P. sp.7 (Indonesia), P. sp (China), Scapliis sp. (Philippine Islands), 
and Peronia sp.2 (Singapore). 
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Discussion 

Morphology. Although several species of Peronia have been described, most information is available only for P. 
verruculata and P. peronii. Nevertheless, P. persiae sp. nov. can be distinguished from all described species by sev¬ 
eral features (see Table SI) and resembles most P. verruculata. Our new species has branchial gills covering only the 
posterior part of the body (Figure 2B) as is described for many Peronia species (Farran 1905; Britton 1984). Only 
P. peronii has gills scattered all over the notum (Plate 1893; Labbe 1934; Solem 1959). Similar to P. verruculata, P. 
persiae sp. nov. lacks spicules in the notum (Awati & Karandikar 1948). Infonnation about this character is miss¬ 
ing for Peronia anomala Labbe, 1934 , Peronia lata Labbe, 1934 and Peronia branchifera Plate, 1893; however, all 
other described Peronia species possess spicules. The visceral cavity in P. persiae sp. nov. is strongly pigmented 
like in P. verruculata (Awati & Karandikar 1948), while the visceral cavity of P. peronii is not pigmented (Labbe 
1934; Plate 1893). Various photoreceptive systems are thoroughly described for P. verruculata (Katagiri etal. 1985) 
including stalk eyes, dorsal eyes, free dermal photoreceptor cells in the notum, and photosensitive neurons in the 
central nervous system. At least the first two types of photoreceptors are present in all Peronia species (Table SI). 
One dorsal tubercle has several eyes and these tubercles are scattered over the dorsal notum. Only P. lata has one 
dorsal eye per tubercle and thus differs from our new species. 

P. persiae sp. nov. lacks the pharyngeal wall teeth, which are present in P. verruculata (Awati & Karandikar 
1948). Additionally, the new species has a larger radula with more rows and more teeth per row, when comparing 
P. verruculata specimens of same size with our specimens (both 65 mm): 71x87.1.87 versus 65><66.1.66 (Awati 


518 ■ Zootaxa 4758 (3) © 2020 Magnolia Press 


MANIEI ETAL. 






& Karandikar 1948). The lateral cusps on both sides of the rachidian teeth are more elongate and pointed, and the 
main cusp is much more developed in P. persiae sp. nov. than in P. verruculata (Plate 1893; Awati & Karandikar 
1948; Britton 1983), thus similar to P. peronii (Labbe 1934). The concave internal part of the lateral teeth was not 
mentioned before (Figure 5C) and might be a special character of the new species. Most onchidiid species have a 
stomach divided into four parts (Weiss & Wagele 1998; Dayrat et al. 2016). Flowever, stomachs with three or two 
parts are also reported (Dayrat 2010a; Dayrat 2010b). Awati & Karandikar (1948) considered the third and fourth 
parts as one single chamber in P. verruculata. The stomach of the here investigated P. persiae sp. nov. has four parts 
in which the fourth has a ciliated epithelium. The intestine with two loops can be assigned to type II. Descriptions 
for P. peronii and P. verruculata vary between the more common intestinal type I (with one loop) and type II (Plate 
1893; Labbe 1934; Marcus & Marcus 1959). 

The position, symmetric or asymmetric, of the nephridium is mentioned as a distinguishing feature for Peronia 
species. P. persiae sp. nov. has an asymmetric nephridium with the right side larger than the left side. According 
to the illustration of P. verruculata (Awati & Karandikar 1948, Figure 38) the nephridium is larger on the left side. 
Flowever, Labbe (1934) mentioned a symmetric nephridium for P. verruculata, as well as for P. peronii. These dif¬ 
ferences in the description show that the nephridium can vary and therefore cannot be considered as a diagnostic 
character. 

In many heterobranch species, a part of the hermaphroditic duct is swollen into an ampulla, in which autosperm 
is stored before release during copulation. An ampulla is not described for any Onchidiidae species except a speci¬ 
men investigated by Marcus (1971) and assigned to Onchidium simrothi Plate, 1893 (Marcus 1971). This specimen 
was recently assigned to the new genus Wallaconchis : Goulding et al. (2018b) re-examined the type material of O. 
simrothi and subsequently considered this name as a nomen dubium. They did not mention an ampulla for any Wal¬ 
laconchis species. It seems that this character is difficult to find in onchidiids, or it is only present in a few species. 
The feature that we describe here in the gonoduct, filled with irregularly lying sperm, is very similar to the structure 
that is described by Schmekel (1971) as ampulla in Nudibranchia. We herein report for the first time the existence 
of ampulla in Peronia. 

The vas deferens often exhibits an enlarged special area (considered as prostate) with glandular cells producing 
prostatic material (e.g., Wagele & Willan 2000). In histological investigations, prostatic cells are visible by their 
typical glandular cells. Interestingly, no prostate was detected in P. persiae sp. nov., while Awati & Karandikar 
(1948) described it in detail for P. verruculata. Unfortunately, no further information about the prostate is available 
for any other Peronia species. The penial complex usually has the same structures in most Onchidiidae species, in¬ 
cluding members of Peronia. The epidermis of the penis is usually covered by spines. We observed some differences 
in P. persiae sp. nov. such as the penial hooks in the anterior and posterior part of the penis being smaller than those 
of the middle region. Short spines in the posterior part are not described in P. verruculata. Additionally, some hooks 
are fork-shaped, and these special forms are scattered all around the penis. The fork-shaped hooks are not known 
for any other Peronia or onchidiid species. Another distinguishable character in P. persiae sp. nov. is the absence 
of calcareous granules which were described on the hooks of the anterior part of the penis in P. verruculata (Awati 
& Karandikar 1948; Britton 1984). 

The penial gland duct of a member of Onchidiidae is investigated histologically here for the first time. Its func¬ 
tion might be provision of substances to help in sperm transfer. No papillae of the penial gland have been observed 
in P. persiae sp. nov., while this structure is described for P. verruculata (Britton 1984) and Peronia gaimardi 
Labbe, 1934 (Labbe 1934). Peronia ferruginea Lesson, 1831 completely lacks the penial gland (Lesson 1830, 
Britton 1984). The colour of the penial needle in P. persiae sp. nov. is light yellow, like in P. verruculata (Awati & 
Karandikar 1948), but dark brown in P. peronii (Plate 1984). The shape of the penial needle opening can vary among 
the Peronia species (Figure 8E). It is shovel-shaped in P. branchifera, oblique or shovel shaped in Peronia gond- 
wana Labbe, 1934, oblique in P. indica, and unilateral thickened in P. verruculata. P. persiae sp. nov. has two new 
different shapes of needle openings, which are not described from any other Peronia species. In most individuals of 
P. persiae sp. nov. it is blunt and only in one individual it is sickle shaped. All other parts of the genital systems are 
similar to the ones described for P. verruculata and P. peronii. 

Based on these combinatory anatomical differences (Table SI), we cannot assign our specimens to any de¬ 
scribed Peronia species and therefore consider it as a distinct new species. 

Molecular data. This first molecular analysis of the Onchidiidae including all available sequences and cover¬ 
ing 11 genera, provide good evidence for the monophyly of Peronia and nearly all other onchidiid genera. Platev- 
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index and Onchidium are not monophyletic, thus supporting former analyses, based on smaller data sets without the 
five recently described genera Wallaconchis, Alionchis, Paromoionchis, Melavonchis and Marmaronchis (Dayrat et 
al. 2017, 2018, 2019; Goulding et al. 2018a; Goulding et al. 2018b). Species delimitation tests clearly show that P. 
persiae sp. nov. is a distinct species. It does not group with any of the other specimens assigned to P. verruculata 
or other, unidentified Peronia sequences. It thus supports the results based on the morphological analysis. Addition¬ 
ally, our results indicate the presence of several further unidentified and potentially undescribed species of Peronia. 
Different species delimitation analyses (GMYC, ABGD, bPTP) based on COI and 16S genes, detected nine species 
within this genus. It is not clear whether these groups represent already described species, because no sequences are 
available for many recognized Peronia species and therefore a revision of the genus is highly needed. 


Conclusions 

Integrative species delimitation analyses are well-established in heterobranchs (e.g. Jorger et al. 2012, Krug et al. 
2013; Padula et al. 2014). Investigating morphology with different methodologies and combining these results with 
those obtained from an array of suitable molecular approaches can provide good evidence in species delimitation 
and descriptions. Our analysis has shown for the first time the presence of a new Peronia species in the intertidal 
flats along the southern coast of Iran. With this new species, ten Peronia species are considered as valid. Future 
analyses will show whether any of the sequences in our analysis confirmed here to belong to the genus Peronia can 
be assigned to one of the former described species. Whether P. persiae sp. nov. is endemic to Iran or has a wider 
geographic distribution is not known, more studies in the Persian Gulf and adjacent areas are necessary. The lack of 
records of other onchidiid species in this region suggests that P. persiae sp. nov. may have evolved certain adapta¬ 
tions to live in the extreme environmental conditions of the intertidal flats of Iran. 
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TABLE SI. Features of valid species within the genus Peronia. Abbreviation in brackets indicates authors and year of publication: A=Awati & Karandikar 1948, At=Attaran etal. 
2015, B=Britton 1984, Bi=Bitaab et al. 2015, Br=Bretnall 1919, C= Cuvier 1804, F=Farran 1905, Fl=FIyman 1999, L=Labbe 1934, Le=Lesson 1830, M=Marcus & Marcus 1959, 
Pa=Patil & Kulkarni 2013, P-PI ate 1893, S=Solem 1959, W=White 1951. 
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FIGURE SI. Rooted ultrametric tree from Bayesian analysis with branches coloured according to the bPTP species delimita¬ 
tion analysis for the genus Peronia based on COI alignment. Bayesian support values are shown above the branches. Blue lines 
indicate branching processes among species, red lines indicate branching processes within species. Wallaconchis graniferus as 
outgroup and identical sequences are removed. 
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FIGURE S2. Rooted ultrametric tree from Bayesian analysis with branches coloured according to the bPTP species delimitation 
analysis for the genus Peronia based on the 16S alignment. Bayesian support values are shown above the branches. Blue lines 
indicate branching processes among species, red lines indicate branching processes within species. Wallaconchis graniferus as 
outgroup and identical sequences are removed. 
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FIGURE S3. Rooted ultrametric tree from Bayesian analysis with branches coloured according to the GMYC species delimita¬ 
tion analysis for the genus Peronia based on COI alignment. Bayesian support values are shown above the branches. Black lines 
indicate branching processes among species, red lines indicate branching processes within species. Wallaconchis graniferus as 
outgroup and identical sequences are removed. 
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FIGURE S4. Rooted ultrametric tree from Bayesian analysis with branches coloured according to the GMYC species delimita¬ 
tion analysis for the genus Peronia based on 16S alignment. Bayesian support values are shown above the branches. Black lines 
indicate branching processes among species, red lines indicate branching processes within species. Wallaconchis graniferus as 
outgroup and identical sequences are removed. 
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